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ABSTRACT 

The angular cross-correlation between two galaxy samples separated in redshift 
is shown to be a useful measure of weak lensing by large-scale structure. Angular 
correlations in faint galaxies arise due to spatial clustering of the galaxies as well 
as gravitational lensing by dark matter along the line-of-sight. The lensing contri- 
bution to the 2-point auto-correlation function is typically small compared to the 
gravitational clustering. However the cross-correlation between two galaxy samples is 
nearly unaffected by gravitational clustering provided their redshift distributions do 
not overlap. The cross-correlation is then induced by magnification bias due to lensing 
by large-scale structure. We compute the expected amplitude of the cross-correlation 
for popular theoretical models of structure formation. For two populations with mean 
redshifts of ~ 0.3 and 1, we find a cross-correlation signal of ~ 1% on arcminute scales 
and ~ 3% on a few arcseconds. The dependence on the cosmological parameters J7 
and A, on the dark matter power spectrum and on the bias factor of the foreground 
galaxy population is explored. 



Key Vifords: galaxies: clustering - cosmology: observations 
large scale structure of the Universe 



gravitational lensing 



1 INTRODUCTION 

The angular auto-correlation function of galaxies has been 
widely used to characterize the large-scale structure in the 
universe (e.g. Peebles 1980). The observed galaxy distribu- 
tion is well described by a power law uj{9) oc d~'' , with slope 
7 ~ 0.8. 

Gravitational lensing by large-scale structure along the 
line-of-sight can alter the amplitude of lu{6) (Gunn 1967). 
Lensing increases the area of a given patch on the sky, thus 
diluting the number density. On the other hand, galaxies 
too faint to be included in a sample of given limiting mag- 
nitude are brightened due to lensing and may therefore be 
included in the sample. The net effect, known as magnifica- 
tion bias, can go either way: it can lead to an enhancement 
or suppression of the observed number density of galaxies, 
depending on the slope of the number-magnitude relation. 
Variations in the number density which are correlated over 
some angular separation alter uo{9). Following Kaiser (1992) 
and Villumsen (1996), Moessner, Jain & Villumsen (1997) 
(henceforth, MJV) have considered the effect of nonlinear 
gravitational evolution and magnification bias on uli{9). MJV 
found that for faint samples with mean redshift 2: ~ 1 lensing 



contributes 5—20% of the signal, depending on the cosmolog- 
ical model and angular scale. Since the lensing contribution 
is small even for distant galaxies, it is diflicult to interpret 
a measurement, especially since it requires knowledge of the 
biasing of high-redshift galaxies relative to the mass. 

In this paper we explore a different statistic, the cross- 
correlation function of two different galaxy samples, in or- 
der to isolate the effect of magnification bias. Consider two 
galaxy samples with non-overlapping redshift distributions. 
If the minimal distance between the two samples is several 
100 Mpc, the effects of gravitational clustering are negligi- 
ble. The cross-correlation function in such a case is affected 
entirely by magnification bias, and the dominant term is pro- 
vided by the lensing effect of dark matter associated with 
the foreground galaxy population. Cross-correlating observ- 
ables affected by gravitational lensing, such as image ellip- 
ticities of galaxies at high redshift, with the positions of 
galaxies at lower redshift has proved to be fruitful for de- 
tecting the effect of gravitational lensing, for example in 
galaxy-galaxy lensing. The cross-correlation of the number 
density of high-redshift quasars with foreground galaxies has 
also been investigated observationally (Benitez & Martinez- 
Gonzalez 1996) and theoretically (Bartelmann 1995, Dolag 
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& Bartelmann 1997). The cross-correlation of magnifica- 
tion and shear has been studied by Kaiser (1992), Sanz, 
Maritnez-Gonzales & Benitez (1997) and Schneider (1997), 
while the auto-correlation of the shear has been computed 
by Blandford et al. (1991), Miralda-Escude (1991), Kaiser 
(1992), Bernardeau, van Waerbeke and Mellier (1996) and 
Jain & Seljak (1997). 

Section 2 provides the formahsm for computing the ef- 
fects of lensing and nonlinear gravitational evolution on the 
angular cross-correlation function of galaxies. Results for 
cold dark matter (CDM)-like models are given in Section 
3. We provide estimates of the errors in observational esti- 
mates of the cross-correlation in Section 4 and conclude in 
Section 5. 



2 ANGULAR CROSS-CORRELATION: 
FORMALISM 

This section provides the formalism for computing the an- 
gular correlation function of galaxies for a given cosmolog- 
ical model and primoridal spectrum of fluctuations for the 
dark matter. The intrinsic gravitational clustering as well as 
the effect of gravitational lensing by the dark matter along 
the line-of-sight are included. Since we consider the angu- 
lar correlation function, only the projected quantities are 
observable. 

We will focus on the cross-correlation of two galaxy 
samples, one at high redshift, z^l, called the background 
sample, and the other at low redshift, z ^ 0.5, called the 
foreground sample. We include in our computations the true 
correlation between galaxies belonging to the high- and the 
low-redshift samples due to the overlap of their redshift dis- 
tributions. As the overlap in the two redshift distributions 
decreases, this intrinsic clustering decreases. The effect we 
wish to isolate is the apparent clustering induced by gravi- 
tational lensing via magnification bias. 

Let ni{(j}) be the number density of galaxies belonging 
to the sample with a low mean redshift (zi), observed in 
the direction (p in the sky, and n2{<j>) that of the sample 
with a higher mean redshift (22) > (.zi)- The angular cross- 
correlation function is then defined as 



u[6) = {5ni{4>)Sn2{i')) 
where 

ni{(t>) — fii 



Sni[cj>) 



(1) 



(2) 



and fii is the average number density of the ith sample. The 
fiuctuation Sui arises due to the true clustering of galaxies 
5nl{(l)), and due to magnification bias Sn^{(l)), 



(3) 



The goal of this paper is to compute the cross-correlation 
io{6) arising from these two efi'ects. 

For simplicity we assume a linear bias model where 
galaxies trace the underlying dark matter fiuctuations. 



5g{x) = b 5{x) 



(4) 



The fluctuations on the sky due to intrinsic clustering are a 
projection of the density fiuctuations along the line-of-sight. 



weighted with the bias factor and the radial distribution 
W{x) of the galaxies 



(5) 



The metric, the comoving radial coordinate x a^nd the co- 
moving angular diameter distance r{x) are introduced in the 
Appendix. We assume for simplicity a constant bias factor 
independent of scale and redshift for each galaxy population. 

To determine the fluctuation due to magnificaton bias 
consider the logarithmic slope s of the number counts of 
galaxies No{m) in a sample with limiting magnitude m (see 
MJV for details) 



dlogiVo(m) 
dm 



(6) 



Magnification by amount fi changes the number counts to 
(e.g. Broadhurst, Taylor & Peacock 1995) 

N'{m)^No{m)fi^-'''-\ (7) 

In the weak lensing limit the magnification is fi = 1 + 2k, 
where the convergence k is a weighted projection of the den- 
sity field along the line-of-sight (see equation ^ below). Since 
K <C 1 for weak lensing, equation ^ for the number counts 
reduces to 



N'{m) = iVo(m) [1 + 5(s - 0.4)^:] 



(8) 



Using g{x) to denote the radial weight function (e.g. Jain & 
Seljak 1997) the convergence k is 



^ Jo °' 



(9) 



The expression for g{x) in terms of the r{x) is given in the 
Appendix. 

Finally, using the above relations we can re-write equa- 
tion ^ for Sni{(f>) as 

5nS)^ r" dxMx)S{r{x)4>,a) , (10) 



with 



Mx) = b^ Wrix) +3nm (2.5s. - 1) 



9^ix) 



(11) 



Inserting the above two equations into Eq. |^, using the small- 
angle approximation ^ <C 1, and assuming that the radial 
weight functions /i(x) vary slowly compared to the scale of 
density perturbations of interest gives (Villumsen 1996), 



u{6) = 47r^ / dxh{x)h{x) 

'0 



X dkkP{x,k)Jo[kr{x)e]S{rix)4>,a) . (12) 
Jo 

The power spectrum of dark matter fiuctuations ^(x, k) is 
defined by 



{5{k)3'{k')) = {27:fP{x,k)5{k-k') 



(13) 



The angular cross- correlation function is composed of 
four terms. In the case of (22) > (zi), these terms are 

+ {5nt{4^5n'i{4>')) + {5nt{4>)5nl{i')) . (14) 
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The first term is due to the intrinsic clustering of the galaxies 
of the two samples where their redshift distributions overlap, 



^ggW = 6i&247r 



dxWi{x)W2ix) 



X / dkkP{x,k)Jo[krix)e] 



(15) 



Ideally we would like this term to be zero, in order to dis- 
tinguish the contribution due to lensing more clearly. This 
could be achieved by obtaining photometric redshifts for the 
galaxies, and selecting two galaxy populations which do not 
overlap in their redshift distributions. The cross-correlation 
of two such samples minimizes the contribution due to in- 
trisic clustering, which removes uncertainties in the predic- 
tions due to unknown physical evolution of the galaxies. 

The second term in equation |lj is due to the lensing of 
the background galaxies by the dark matter in front of it, 
which is traced by the foreground galaxies. The correlation 
thus induced between galaxies in the two samples is given 
by, 



bimm{2.5s2 - l)47r^ 



X / dkkP{x,k)Jo[kr{x)0] 
lo 



(16) 



The third term is due to dark matter in front of both of 
the galaxy samples doing the lensing. The fourth is due to 
dark matter traced by the background galaxies lensing the 
foreground galaxies. It is non-zero only if there is an overlap 
in the redshift distributions of the two samples. For all cases 
of interest these two terms are negligible compared to the 
second term, ujgi. 

2.1 Dependence on the cosmological model 

Equation 15 shows how uj{9) depends linearly on flm, aside 
from the dependences on f2,„ and Qa contained in the line- 
of-sight integral. These arise from two sources: (i) the dis- 
tance factors contained in Jo, g{x) ^i^d W{x), and (ii) the 
growth and amplitude of the power spectrum. In the linear 
regime, P{Xj k) depends on the linear growing mode of den- 
sity perturbations -D(x) the normalization ag (which in 
turn can depend on the cosmological parameters) as 

P{x,k)^[asD{x)f . (17) 

The linear growing mode is well approximated by (Carroll, 
Press & Turner 1992) 

5 

— I 

2 

+ (1 



D{X) 



■afl{a) [n{a)*^'' - \{a) 

-fi(a)/2)(l + A(a)/70)] 



(18) 



where we have defined, following Mo, Jing & Borner (1996), 
the time dependent fractions of density in matter and vac- 
uum energy, 57(a) and A(a) in terms of the present-day values 
f2m and Oa, 



and 



A(a) 



, + S7m(l — a) + SlA(a3 — a) 



a + f2m(l — a) + flA{a^ — a) 



(19) 



(20) 



Moreover, the spatial geometries differ in different models, 
leading to a dependence of the angular distance r{x) on $1™ 
and flA according to Eq. ^ 

The redshift distribution of galaxies can be modelled by 



^gr[3//3] 



exp [-{z/zqY 



(21) 



for P = 2.5, which agrees reasonably well with the redshift 
distribution estimated for the Hubble deep field from photo- 
metric redshifts (Mobasher et al. 1996). The mean redshift 
is then given by 

The four different cosmological models we consider are 
a fiat universe with Q,m = 1, an open model with Q.m. = 
0.3, and a flat A— dominated model, with Q.m = 0.3 and 
JIa = 0.7. These three models are normalized to cluster 
abundances (White, Efstathiou & Frenk 1993; Viana & Lid- 
die 1995; Ecke, Cole & Frenk 1996; Pen 1996), 



0-8 



0.6 r2: 



(23) 



(For r2,„ = 0.3, we use erg — 1, which is close to the results of 
the approximate formula given above.) Our fourth model is 
a flat universe, Qm = 1, with a high normalization of erg = 1. 

The nonlinear power spectrum is obtained from the lin- 
ear one through the fitting formulae of Jain, Mo & White 

(1995) for Q,rn = 1, and from those of Peacock & Dodds 

(1996) for the open and A— dominated models. These fit- 
ting formulae are based on the idea of relating the nonlinear 
power spectrum at scale k to the linear power spectrum at 
a larger scale kL, where the relation k{kL) depends on the 
power spectrum itself (Hamilton et al. 1991). They have 
been calibrated from and tested extensively against N-body 
simulations. For the linear power spectrum we take a CDM- 
like spectrum with shape parameter F = 0.25, which pro- 
vides a good fit to observations. The transfer function for 
the initially scale-invariant power spectrum is taken from 
Bardeen et al. (1992). 



3 CROSS-CORRELATION FUNCTION FOR 
CDM-LIKE POWER SPECTRA 

For samples not overlapping in their redshift distributions, 
the cross-correlation cogi induced by lensing dominates the 
contribution to the cross-correlation function. The signal 
due to intrinsic clustering is negligibly small, though if 
the redshift distributions overlap significantly, it would 
swamp the lensing signal. In this section we present results 
on the cross-correlation of two galaxy samples with non- 
overlapping redshift distributions. 

We computed ujgi{0 — 1') for two populations with 
mean redshifts of 0.3 and 1.5 using the redshift distribution 
n{z) of Eq. We found that the results differ by 5 - 10% 
from those obtained using a delta-function redshift distri- 
bution. Therefore we will use the simpler form of a delta- 
function distribution for computing u!gi{9). For calculating 
contributions due to intrinsic clustering, however, it is nec- 
essary to use the full n{z) of Eq. We do this in Section ^ 
below to estimate the relative error made by misidentifying 
background galaxies as foreground ones. 
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Figure 1. The cross-correlation function cjjj as a function of 6 is shown for the four cosmological models. Since the number count slope 
s = 0.2, magnification bias induces an anti-correlation between the foreground and background sample. Hence the sign of LOgi is negative. 
The differences between the four models are discussed in the text. 



In Figure |l| we plot Ugi as a function of 6 for a fore- 
ground galaxy sample at zi = 0.3 and a background sample 
at Z2 = 1.5, for the four cosmological models described in 
section 2.1. We choose a bias factor of & = l/cg, which is 
in agreement with large-scale galaxy clustering data; Ugi 
is proportional to b. We assume a number count slope of 
a = 0.2 for the background sample. For s < 0.4 the induced 
correlations are negative since uigiUgi oc (s — 0.4). This linear 
relation also makes it simple to scale our results to other val- 
ues of s. A sample with a slope close to 0.2 may be obtained 
by defining color selected subsamples, using the fact that 
the number count slope is a decreasing function of — / 
color (Villumsen, Freudling & da Costa 1996; Broadhurst 
et al. 1997). The price of selecting a subsample is a smaller 
number of galaxies and therefore larger Poisson errors; so 
in practice a careful cut suited to the available data would 
need to be made. 



Figure |l| shows that the typical cross-correlation signal 
expected on sub-arcminute scales is a few percent. On an- 
gular scales larger than 1 arcminute the signal drops to less 
than a percent. For fixed ag it is largest for the Einstein- 
de Sitter model and smallest for the cosmological constant 
model, at least on arcminute scales or smaller. However if ug 
is determined from cluster-abundances, the cross-correlation 
on small scales is largest for the open model. 

The dependence of the cross-correlation on the redshift 
of the background sample is shown in Figure ^. We have 
plotted the cross-correlation function uogiifl = 0.1') for a 
foreground galaxy sample at zi — 0.3 as a function of the 
redshift Z2 of the background sample. Figure ^ is the same 
as Figure ^, but for 9 = 0.2'. These figures show the slow 
increase in the amplitude of the signal with (2:2) above a 
redshift of 1. There is no significant variation in the shape 
of the curves among the four cosmological models. Thus if 
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Figure 2. Cross-correlation function as a function of Z2 is shown 
for zi = 0.3 at 6 = 0.1' for the four cosmological models. The neg- 
ative correlations induced by magnification bias become stronger 
with Z2, but do not change much beyond a redshift of 1. 



the amplitudes were normalized to the same value at 2:2 — 1, 
there would be very little difference between the curves at 
higher Z2. 

The difference between the predictions of the four mod- 
els shown in Figures can be qualitatively understood as 
follows. The dominant dependence arises due to the factor 
of Qm outside the integral in equation |l^ for uigi. Taking 
the normalization of the power spectrum into account this 
is reduced to O^'* for the cluster-abundance normalization 
of erg and the bias relation b — l/ag. This is because the 
expression for ujgi depends explicitly on the factor bflmi^s- 

The line-of-sight integral in equation ^ further weakens 
the dependence on Qrn- In a low-flm universe, the growth of 
perturbations is slowed down at late times. Hence, normaliz- 
ing to present-day cluster abundances leads to a higher nor- 
malization at earlier times compared to the ^Im = 1 models. 
This in turn means that nonlinear effects, which are signifi- 
cant on angular scales of 1' or less, become important earlier 
on and lead to a larger enhancement due to nonlinear clus- 
tering by today. The nonlinear enhancement is reinforced 
by a geometrical effect: for lower Qm, and even more so for 
larger 57a (for given Qm), the physical distance to a given 
redshift is larger. This leads to a larger lensing path-length 
and thus a further increase in the lensing signal. 

The combination of all these effects is shown in Fig- 
ure |l[ On scales well below an arcminute the signal for the 
open model becomes comparable to that in the Einstein-de 
Sitter model. This is due to the dominance of the nonlinear 
enhancement on these scales. The curve for the open model 



Figure 3. Cross-correlation function as a function of 22 is shown 
for 6 = 0.2' for the four cosmological models. All parameters 
except 9 are as in the previous figure. 



is also distinctly steeper than for the others. For 9 > 2' , 
however, there is no significant variation in the shape of the 
curves. For the A— dominated model the effect of the growing 
mode is not as strong as for the open model. Thus on small 
scales, even though the geometric effect gives a stronger en- 
hancement than for the open model, the net amplitude of 
uigi is smaller than for the open model. 

Finally, note that cUgi is proportional to the bias factor of 
the foreground sample. If this bias factor is larger or smaller 
than the value of 6 = l/os which we assumed to fit the large- 
scale structure data, Ugi will be correspondingly altered. 



4 ESTIMATE OF ERRORS 

We consider two sources of error involved in an observa- 
tional determination of the cross-correlation from two galaxy 
samples. The first potential error arises when a background 
galaxy is mis-identified as a foreground galaxy. In this case 
the auto-correlation of the background galaxy sample will 
erroneously enter into the measured cross-correlation. Say 
e% of the background galaxies' redshifts are sufficiently mis- 
estimated that they are taken to be part of the foreground 
sample. The observed cross-correlation is then given by 

N2 

IjJm ^ UJgl + ';j^UJgg{z2) , (24) 

where the second term is the error made due to assigning 
galaxies to the wrong sample. It is proportional to the frac- 
tion of galaxies which is mis-identified and to ujgg{z2), the 
auto-correlation function evaluated at the mean redshift 
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Figure 4. True cross-correlation function Ugi (lower set of matching curves) and measured cross-correlation function u)m (upper set of 
matching curves) when e% of the background galaxies are misidentified as foreground ones, as a function of 22, for 9 = 0.1', 21 = 0.3 
and the four cosmological models. 



of the background sample, for the redshift distribution of 
Eq. ^ We have calculated the auto-correlation at 1' and 
scaled it to 0.1' assuming a power law slope of —0.8; TVi 
and N2 are the number of galaxies in the foreground and 
background samples. 

In Figure ^ we plot the cross-correlation Ljgi{6 = 0.1') 
(lower set of matching curves) together with the mea- 
sured cross-correlation function Um (upper set of matching 
curves). Two values of e, 1% and 5% are used (for samples 
of equal size). The results show that for 22 > 1 the error is 
small for e = 1% but not when e = 5%. This sets an approx- 
imate standard required for using photometric redshifts or 
other possible methods to select the two galaxy samples. 

A second source of error is the statistical uncertainty 
in estimating the angular correlations. Using a Poisson dis- 
tribution to estimate the error in u) provides a rough guide 
for the number of galaxies required to estimate the cross- 
correlation signal. The standard deviation 5io{9) in the es- 
timate of uj{9) for a random distribution of objects is given 
by (Peebles 1980) 



5Lo{ef = ^ 



N1N2 sn 



(25) 



where fl is the solid angle subtended by the survey area, and 
5Q is the fraction in the bin used for angle 9. Note that dcj^ 
is just the inverse of the number of pairs in a given bin in 9. 

For a sample with about 10^ galaxies per 0.01 degree"^ 
(e.g. as in each field of Woods & Fahlman (1997) whose 
sample reaches limiting magnitudes of V ^ 25, R ~ 25, 
7 ~ 24), the above estimate gives iJaj ~ 4 x 10""^ if 10 bins in 
9 are used. Thus in excess of about 10^ galaxies each in the 
foreground and background sample would be required for 
a detection of a ^1% cross-correlation signal with a high 
level of significance. 



5 CONCLUSIONS 

We have presented results for the cross-correlation of two 
galaxy samples with different redshift distributions. The sig- 
nal is dominated by the effect of magnification bias due to 
weak lensing. To ensure that the contribution from gravi- 
tational clustering of the galaxies is negligible, we have as- 
sumed that the redshift distributions of the two samples do 
not overlap. With the use of photometric redshifts (e.g. Con- 
nolly et al. 1995; Sawicki, Lin & Yee 1997), it is feasible to 
obtain deep galaxy samples that can be separated into sub- 
samples with desired redshift distributions. If only limiting 
magnitudes are used to create two sub-samples, then there 
will be a significant overlap and the interpretation of the 
signal is not as clean. Theoretical predictions can however 
be made for the expected signal assuming a redshift distri- 
bution. 

The results shown in Figures 1-3 demonstrate that most 
models predict a signal of 1-4% for the cross-correlation 
function. These numbers apply for a background sample 
with a mean redshift ^ 1 and a number count slope of 0.2, 
on angular scales from a few arcseconds to an arcminute. 
As argued in section 4, the measurement of such a signal 
appears feasible in the near future. 

The cross-correlation function is a measure of the pro- 
jected dark matter power spectrum. For a given spectrum, 
the variation with angle on small scales is largest for open 
cosmological models and thus provides a probe of Q. The 
cross-correlation however is also proportional to the bias fac- 
tor of the foreground {z ~ 0.3 — 0.5) galaxy sample. Thus 
there is a degeneracy in the dependance on the cosmologi- 
cal model and the biasing of intermediate redshift galaxies. 
The bias factor can vary with the redshift of the foreground 
sample and with SI. Given a model or empirical measure 
of the bias factor, the cross-correlation can be used as a 
probe of the power spectrum and Q. Else, for a given cos- 
mology, it can constrain the bias factor of galaxies. Ideally, 
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by combining the cross-correlation with other lensing mea- 
sures which are independent of bias, such as the ellipticity 
auto-correlation function, constraints on the cosmological 
model as well as on the biasing of galaxies at intermediate 
redshifts can be obtained. 
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APPENDIX A: NOTATION 

Following the notation of Jain & Seljak (1997) we introduce 
the unperturbed metric 



da 



\t) {-dr'^ + dx^ + r^(d6l^ + sin^ 9d<fP)) , 



(26) 



with r being conformal time, and x the radial comoving 
distance, xh is used to denote the distance to the horizon. 
The comoving angular diameter distance r{x) is, 

{K-^^'^sinK^/^X, K >0 
X, K = (27) 
{-K)-^^^smh{-K)'-^^X, K<0 

where K is the spatial curvature given hy K — -~Ho{l — 
flrn — Q — A) with Ho being the Hubble parameter today. 

With VF(x) denoting the radial distribution of galaxies 
in the sample, the radial weight function g{x) is given by 



9ix) = r{x) '^^W{^)d^ . (28) 

For a delta-function distribution of galaxies, W{x ) = 
5d[x' - Xs), and g[x) reduces to g[x) = rixVixs - 
x)/r{xs)- 
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